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ABSTRACT 

Galaxy cluster Abell 3827 hosts the stellar remnants of four almost equally bright 
elliptical galaxies within a core of radius 10 kpc. Such corrugation of the stellar dis¬ 
tribution is very rare, and suggests recent formation by several simultaneous mergers. 
We map the distribution of associated dark matter, using new Hubble Space Telescope 
imaging and VLT/MUSE integral field spectroscopy of a gravitationally lensed system 
threaded through the cluster core. We find that each of the central galaxies retains a 
dark matter halo, but that (at least) one of these is spatially offset from its stars. The 
best-constrained offset is 1.621 q 49 kpc, where the 68% confidence limit includes both 
statistical error and systematic biases in mass modelling. Such offsets are not seen 
in field galaxies, but are predicted during the long infall to a cluster, if dark matter 
self-interactions generate an extra drag force. With such a small physical separation, it 
is difficult to definitively rule out astrophysical effects operating exclusively in dense 
cluster core environments - but if interpreted solely as evidence for self-interacting 
dark matter, this offset implies a cross-section (Tdm/hi ^ (1.7 ± 0.7) x 10“^cm^/g 
X (tinfaii/10^ yrs)“^, where tinfaii is the infall duration. 

Key words: dark matter — astroparticle physics — galaxies: clusters: individual: 
Abell 3827 — gravitational leasing: strong 
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1 INTRODUCTION 

Many lines of evidence now agree that most mass in the 
Universe is in the form of dark matter, which interacts 
mainly via the force of gravity. The identity and detailed 
phenomenology of dark matter remain poorly understood. 
However, its gravitational attraction pulls low-mass systems 
into a series of hierarchical mergers through which dark and 
ordinary matter are gradually assembled into giant clusters 
of galaxies ( Davis et al.|1985 ). 

The typically smooth distribution of light in galaxy 
clusters visible today shows that merging systems have their 
gas content efficiently removed into the intra-cluster medium 
by ram pressure stripping, even while they pass the virial 
radius (Smith et al. 2010| Wu et al. 2012). The longevity 
of accompanying dark matter is less well understood, but 
the timescale for its dissipation is a key ingredient in semi- 
analytic models of structure formation ( [Dariush et al.|2010 ). 
Full numerical simulations predict that the dark matter is 


eventually smoothed (Gao et al. 2004 Nagai et al. 2005 
Bahe et al.|2012 ), but disagree about the timescale and the 


radius/orbits on which stripping occurs (Diemand, Kuhlen 
Madau||2007| [Peharrubia, McConnachie k. Navarro||2008| 


Wetzel, Cohn & White 2009). Observations have shown that. 


as L* galaxies enter a galaxy cluster from the field, tidal 


gravitational stripping of their dark matter (Mandelbaum 


et al. 


2006 


2007 


Parker et al. 
their masses by ~ 10^"^ M© to 


Gillis et al. 2013) reduces 
from a radius of 


10 " 


5Mpc to 1 Mpc (Limousin et al. 2007 2012 Natarajan et 


al. 2009). This stripping occurs at a rate consistent with 


simulations, but has not been followed to the central tens 
of kiloparsecs, which is where the predictions of simulations 
disagree. 

Mergers of dark matter substructures into a massive 
galaxy cluster also reveal the fundamental properties of dark 
matter particles. The different non-gravitational forces act¬ 
ing on dark matter and standard model particles have been 
highlighted most visibly in collisions like the ‘Bullet Clus- 


ter’ 1E0657-56 ( 

Glowe et al.|2004 

[2006| 

Bradac et al. 

2006), 

Abell 520 (Mahdavi et al. 

2007 1 

1 Glowe et al. 

2012 

Jee et 

ZlpONl, MACSJ0025-12 ( 

Bradac et al.[[2008 

1, Abell 2744 


al.||2012 ). Infalling gas (of standard model particles) is sub¬ 


ject to ram pressure and tends to lag behind non-interacting 
dark matter ( Young et al.||20lT ). Measurements of this lag 


yielded an upper limit on dark matter’s self-interaction cross 
section cr/m < 1.2 cm^/g if the particle momentum exchange 


is isotropic, or cr/m<0.7 cm^/g if it is directional (Randall 
et al.|200^ [Kahlhoefer et al.|2014| . More interestingly still. 


if dark matter has (even a small) non-zero self-interaction 
cross-section, infalling dark matter will eventually lag be¬ 
hind old stars ( Massey, Kitching fc Nagai|2011[ Williams & 


Saha|2QlT] Kahlhoefer et al.|2014 Harvey et al.|2013 2014). 
Self-interactions confined within the dark sector can poten¬ 
tially have much larger cross-sections than those between 
dark matter and standard model particles, which are con¬ 
strained by collider and direct detection experiments {e.g. 
Peter et al.||2012 ). 

The galaxy cluster Abell 3827 (RA=22h01'4971, 
Dec=—59° 57'15", z=0.099. X-ray l uminosity Lx ^ ^ ^ 
10^^erg/s in the 0.1-2.4keV band, De Plaa et al. 2007) 


it hosts the remnant stellar nuclei of four bright elliptical 
galaxies within the central 10 kpc. Such corrugation of the 


stellar distribution is very rare: only Abell 2261 (Goe et al. 


2012 


2012 


Postman et al. 2012) and MAGSJ0717 (Limousin et al. 
Jauzac et al.||2012 ) are even comparably corrugated. 


All three clusters are still forming, through several simulta¬ 
neous mergers - and can be used to investigate the late-stage 
dissipation of dark matter infalling through the same envi¬ 
ronment. Moreover, Abell 3827 has a unique strong gravi¬ 
tational lens system threaded between its multiple central 
galaxies ( Garrasco et al.|2010| . This enables the distribution 
of its otherwise invisible dark matter to be mapped (for re¬ 
views of gravitational leasing, see [Ba£tehnnmj j^Sctm^^ 


2001||Refregier|2003l |Hoekstra &: Jain| 2008| [Massey, Kitch¬ 

ing Richard||2010| |Bartelmann||2010] [Kneib &: Natarajan 
2011| ). The cluster even lies within the optimum redshift 
range 0.05 < z <0.1 to measure small physical separations 


between dark and ordinary matter ( [Massey, Kitching Na- 
gai||2011| 


Indeed, ground-based imaging ( Williams Sah^|2011 


Mohammed et al. 2014) suggests that dark matter asso¬ 


ciated with one of the central galaxies in Abell 3827 (the 
one where its position is best constrained, ‘nucleus’ N.l) is 
offset by ~ 3" (6 kpc) from the stars. Such offsets are not 
seen in isolated field galaxies {e.g. Koopmans et ah] 2006| 
[Gavazzi et al.|[2007| ). Interpreting the offset via a model in 
which tinfaii is the time since infall, implies a lower limit of 
cr/m >4.5 X 10“®(tinfaii/10^°yr)“^ cm^/g. This is potentially 
the first detection of non-gravitational forces acting on dark 
matter. 

In this paper, we present new Hubble Space Telescope 
(HST) imaging and Very Large Telescope (VLT) integral- 
field spectroscopy to hone measurements of the dark matter 
distribution. We describe the new data in section and 
our mapping of visible light/dark matter in section^ We 
describe our results in section and discuss their impli¬ 
cations in section We conclude in section Through¬ 
out this paper, we quote magnitudes in the AB system and 
adopt a cosmological model with Dm = 0.3, Da = 0.7 and 
Ho = 70km/s/Mpc, in which l" corresponds to 1.828 kpc 
at the redshift of the cluster. 


2 DATA 

2.1 HST imaging 

We imaged the galaxy cluster Abell 3827 using the Hub¬ 
ble Space Telescope (HST), programme GO-12817. Obser¬ 
vations with the Advanced Camera for Surveys (ACS)/Wide 
Field Channel (WFC) during October 2013 consisted of 
5244 s in optical band F814W (in the core, with half 
that depth across a wider area) and 2452 s in optical 
band F606W. Observations with the Wide Field Camera 

3 (WFC3) in August 2013 consisted of 5871s in UV band 
F336W and 2212 s in near-IR band F160W. 

The raw data exhibited spurious trailing due to Gharge 
Transfer Inefficiency in imaging detectors damaged by radia¬ 
tion. We corrected this trailing using the so ftware ofjMassey 


et al. (2010 


din (2010) 


2014) for the ACS/WFC and [Anderson B^ 


Anderson (2014) for WFC3/UVIS. Subsequent 


is particularly interesting for substructure studies because data reduction then followed the standard procedures of 
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Figure 1. Hubble Space Telescope image of galaxy cluster Abell 3827, showing the F160W (red), F606W (green) and F336W (blue) 
bands. The colour scale is logarithmic. Labels show the four bright (plus one faint) central galaxies, foreground stars and background 
lensed galaxies. An object previously referred to as N.5 is actually a star. All these identifications are confirmed spectroscopically. 


CALACS v2012.2 ([Smith et al.|2012|| and calwfS v2.7 (|Sabbi 


et al.|2009|. We stacked individual exposures using DRIZZLE 


Fruchter fc Hook|2002 ) with a Gaussian convolution kernel 
and parameter pixfrac=0.8, then aligned the different ob¬ 
servations into the common coordinate system of the F814W 
data using tweakback. Figureshows a multicolour image 
of the cluster core. 


2.2 VLT spectroscopy 

We first obtained spectroscopy across the cluster core using 
the VLT/ VIMOS integrated field unit (IFU; |Le Fevre et al. 


2003 2013), programme 093.A-0237. Total exposure times 


were 6 hours in the HR-blue filter (spanning a wavelength 
range 370-535 nm with spectral resolution A/AA = 200) dur¬ 
ing July 2014 and 5 hours in the MR-orange filter (490- 
1015 nm, with A/AA=1100) during August 2014. All obser¬ 
vations were obtained in photometric conditions and < 0(^6 
seeing, using the 27"x 21" field of view; in this configuration, 
each pixel is 0('66 on a side. 


Since the cluster core is high surface brightness across 
the entire VIMOS field of view, we interspersed every three 
exposures on target with one offset by ~ 2' to record (and 
subtract) the sky background. The three on-source expo- 
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sures were dithered by 1^/3 (2 pixels) to account for bad fibres 
and cosmetics. 

To reduce the raw data we used the VIMOS esorex 
pipeline, which extracts the fibres, wavelength calibrates and 
Hatfields the data, and forms the data cube. We used the 
temporally adjacent sky exposure to perform sky subtrac¬ 
tion (on a quadrant-by-quadrant basis), then mosaiced all 
of the exposures using a clipped average (using the bright 
stars to measure the relative offsets between cubes). We con¬ 
structed (wavelength collapsed) continuum images from the 
cubes and aligned the cube to the HST imaging, then ex¬ 
tracted spectra for each of the continuum sources. To search 
for emission from the strong lensing features, we applied a 
mask to the cube and extracted both one- and two- dimen¬ 
sional spectra. Due to the different resolutions of the HR- 
blue and MR-orange observations, we analysed the final two 
data cubes separately, but overall they provide a continuous 
wavelength coverage from 370 to 1015 nm. This is perfectly 
sufficient for our analysis of the cluster light. 

The lensed galaxy threaded through the cluster core 
(labelled A in figure was originally identified in Gemini 
imaging by Carrasco et al. (2010), who also used long-slit 
Gemini spectroscopy to obtain a redshift z = 0.204. However, 
our IFU spectroscopy did not confirm this. We instead found 
only one bright emission line at 835.5 nm, which is not asso¬ 
ciated with lines at the foreground cluster redshift and whose 
2D morphology traces the lensed image. The emission line 
could have been Ha at z ^ 0.27 or [Oil] at z ^ 1.24 — but 
the low resolution and signal-to-noise of the VIMOS data 
precluded robust identification (the lack of other features, 
such as [Nil] A6853 for ; 2 ; = 0.27 or [Oiii] A5007 for 2 : = 1.24 
may have been due to low metallicity and / or low signal- 
to-noise in the lines). Moreover, the flat-field (in)stability 
in VIMOS data left strong residuals after subtracting fore¬ 
ground emission from the cluster galaxies. In particular, this 
made it difficult to robustly determine the morphology of the 
arc near N.l, and — as we shall see in section [T^ — this 
provides the most diagnostic power in the lens modelling. 

To confirm the redshift of lensed system A, and to mea¬ 
sure its morphology around N.l, we observed the cluster core 
with the VLT Multi-Unit Spectroscopic Explorer (MUSE) 
IFU spectrograph ( Bacon et al.|201(3 ) during Director’s Dis¬ 
cretionary Time in December 2014, program 294.A-5014. 
The awarded exposure time of 1 hour was split in to 3 x 1200 s 
exposures, which were dithered by ~10'' to account for cos¬ 
mic rays and defects. These observations were taken in dark 
time, <0.7" U-band seeing and good atmospheric trans¬ 
parency. We stacked these with an extra 1200 s exposure that 
was taken in good seeing during twilight, which marginally 
improves the signal to noise. MUSE has a larger 1' x 1' field 
of view and excellent Hatfield stability, so no extra sky ex¬ 
posures were required. The data were reduced using vl.O of 
the ESOREX pipeline, which extracts the spectra, wavelength 
calibrates, HatHelds the data and forms the data cube. These 
were registered and stacked using the exp_COMBINE routine. 
The (much) higher throughput and higher spectra resolution 
(A/AA = 3000) of MUSE yielded greatly improved signal-to- 
noise in both continuum and emission lines (see Appendix [A| 
for a comparison). For all our analysis of the background 
sources, we therefore use only the MUSE data. 

In lensed system A, the MUSE data resolved the 
835.5 nm emission line as the [Oil] A3726.8, 3729.2 doublet at 
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Figure 2. Observed spectra at the locations of background galax¬ 
ies A (top panel) and B (bottom panel), smoothed for clarity with 
a Gaussian of width 5 A. In both cases, the grey line shows the 
spectrum of nearby emission from the foreground cluster, spa¬ 
tially interpolated to the position of the background galaxy. The 
coloured insert shows galaxy B’s 2D velocity field in a 6" x 6" 
region, for all IFU pixels where Ha emission is detected at signal 
to noise > 4. It indicates a rotationally supported disc. 



redshift za = 1.24145 ± 0.00002, confirmed by the additional 
identification of Mgll absorption at 627.0 nm (figure]^. The 
MUSE data also allowed a much improved continuum image 
to be constructed around the [Oil] emission, then subtracted 
to leave a higher fidelity [Oil] narrow-band image (figure [3|. 

For arc B, the MUSE data confirms [Carrasco et al.| 
( 2010| ’s Gemini long-slit redshift, finding zb = 0.4082 ± 
0.0001 at the centre of the arc. Two blue knots at the north 
end dominate the [Oil] and [Olll] emission, but Ha emis¬ 
sion is visible across its entire length (all three features are 
at the same redshift). A 2D map of the best-fit wavelength 
of the Ha emission shows a monotonic velocity gradient of 
~ 200km/s from north to south (figure]^. We are confident 
that no multiple images are present, having inspected both 
the HST imaging and narrow-band images created from the 
MUSE data at the wavelengths of the emission lines. 
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Figure 3. Hubble Space Telescope image of the core of Abell 3827, with light from the four bright galaxies subtracted to reveal the 
background lens system. Colours show the F814W (red), F606W (green) and F336W (blue) bands, and the colour scale is square root. 
Linearly spaced contours show line emission at 835.5 nm from VLT/MUSE integral field spectroscopy. Residual emission near galaxies 
N.3 and N.4 may be a demagnified fifth image or merely imperfect foreground subtraction, so we do not use it in our analysis. 


3 ANALYSIS 

3.1 Modelling the cluster light distribution 

It is apparent from the high resolution imaging (figure]^ and 
our IFU spectroscopy that Abell 3827 contains four bright 
central galaxies, N.1-N.4. The object labelled N.5 by |Car-| 
rasco et al. (2010) is a Milky Way star: it is a point source 


in the HST imaging, and its spectrum contains z = 0 Call 
H and K absorption lines that are not present in adjacent 
sources. Their spectroscopy of N.5 was probably contami¬ 
nated by the nearby bright cluster galaxies and diffuse intra¬ 


cluster light. On the other hand, the westernmost object that 


Carrasco et al. (2010) identihed as a star, is actually a faint 


cluster member galaxy at z = 0.1000 ± 0.0002. To avoid 
confusion, we denote this galaxy N.6. 


In the optical HST imaging, we use galfit (Peng et 


al. 2010) to simultaneously ht the light distribution from 


the four bright galaxies and the two stars. Most galaxies are 
weh-ht by a single component model with a Sersic prohle, 
although a double-component model using two Sersic pro- 
hles (with the same centre) is preferred for galaxy N.3 (and 
galaxy N.2 in the F814W band). Positions in the F814W 
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Table 1. Total integrated flux of the bright central galaxies, and their derived stellar masses. In both HST/ACS bands (independently), 
we use GALFIT ( |Peng et al.|2010| to simultaneously fit the emission from all four galaxies. Asterisks denote that a two-component (co¬ 
centered) Sersic model was preferred. Positions are listed from the F814W band, and are consistent with those from the F606W band. 
Stellar masses M* interpret the single-band AB magnitude flux via the models of |Bruznal &: Chariot] ( [?003| >, assuming a |Chabri^ ( |200^ 
IMF, solar metallicity, and formation redshift Zf = 3. Redshifts, 3cr upper limits on Ho; flux (interpreted as limits on star formation rate 
via |Kennicntt|I998[ also converted to a |Chabrier|2003| IMF), and stellar velocity dispersions are measured from VLT spectroscopy. 



RA 

Dec 

z 

F814W band 
mag M* [Mq] 

F606W band 
mag M* [Mq] 

Ha flux 
[erg/s/cm^] 

SFR 

We/yv] 

[km / s] 

N.l 

330.47518 

-59.945997 

0.09891T0.00032 

16.86 

1.04x10“ 

17.52 

1.00x10“ 

< 6.04x10-1® 

< 0.14 

332 

N.2 

330.47233 

-59.945439 

0.09928T0.00017 

15.74* 

2.92x10“ 

16.55 

2.46x10“ 

< 2.59x10-1® 

< 0.06 

377 

N.3 

330.46978 

-59.944903 

0.09973T0.00016 

15.69* 

3.06x10“ 

16.42* 

2.77x10“ 

< 6.31x10-1® 

< 0.14 

326 

N.4 

330.46999 

-59.946322 

0.09636T0.00026 

16.18 

1.94x10“ 

16.73 

2.08x10“ 

< 1.57x10-15 

< 0.26 

192 


band are listed in table those in F606W are consistent 
within O'.'OOG for N.1-3 and, OVO 6 I for N.4, due to its prox¬ 
imity to a bright star. To model emission from the stars, 
we shift and rescale an isolated star in the same image. The 
best-fit galaxy fluxes do not depend significantly upon the 
choice of isolated star or the alternative use of a TinyTim 
model star ( Krist, Hook Stoehr|20li] Rhodes et al.|2007 ). 
The photometric errors are dominated by our assumption 
of analytic functions to fit the light profiles. The fluxes are 
likely to be an upper limit because they are computed by 
integrating these analytic functions to inhnite radius, and 
may also include a component of diffuse intra-cluster light. 

In the VLT/VIMOS spectroscopy, we measure the red- 
shift of galaxies N.1-N.4 and N.6 by htting a Gaussian to 
the H, K and G-band absorption features. None of the clus¬ 
ter galaxies exhibits Ha line emission, although we attempt 
to ht a Gaussian at its redshifted wavelength to obtain 3cr 
upper limits on the Ha flux. To measure the stellar veloc¬ 
ity dispersion, we cross-correlate our spectra with broadened 
stellar templates from Vazedkis (1999). These measurements 
are presented in table^ 


3.2 Strong lens identifications 

Figure presents a multicolour image of Abell 3827, after 
subtracting the best-ht model of optical emission from the 
optical HST bands to reveal the morphology of the gravita¬ 
tionally lensed system. In the UV HST imaging, the contrast 
between cluster member galaxies and the background lens 
system is much lower, so we do not ht and subtract the 
foreground flux. 

Lensed image A is an almost face-on spiral galaxy, with 
a bulge and many resolved knots of star formation that can 
ah be used as independent lensed sources. The association 
of knots between multiple images is not perfectly clear, due 
to the bright intra-cluster light and a surprising density of 
point sources, particularly near galaxy N.l; we present the 
most likely identihcations in table but analyse alternative 
conhgurations in appendix ^ Gontours in hgure show a 
(continuum subtracted) narrow-band image created by sub¬ 
tracting continuum emission from each spatial pixel (ht- 
ted using a low-order polynomial over the wavelength range 
830-840 nm) then collapsing the MUSE IFU data cube over 
zbSOOkms”^ from the peak of the emission. The 2D map 
of this line emission matches precisely the lensed galaxy’s 
broad-band morphology. Variations in the intensity of the 


Table 2. Locations of multiply imaged systems. Images Ao.n are 
the bulge, and images A[a-f].n are knots of star formation in the 
spiral arms. Index n is sorted in order of arrival time according 
to our fiducial model (see table [^. Columns denote the ID and 
position of the image, its major and minor axes, and the angle of 
its major axis on the sky, anticlockwise from west. 


Name 

RA 

Dec 

Major 

Minor 

Angle 

Ao.l 

330.47479 

-59.943580 

0"33 

0722 

25° 

A 0.2 

330.46649 

-59.946650 

(X'35 

0723 

75° 

A 0.3 

330.46828 

-59.944112 

0"43 

0716 

140° 

A 0.4 

330.47407 

-59.946239 

(X'39 

(X'25 

85° 

Aa.l 

330.47559 

-59.944009 

0"16 

0714 

151° 

Aa.2 

330.46725 

-59.947321 

(X'16 

(X'14 

140° 

Aa.3 

330.46871 

-59.944215 

0"16 

0714 

131° 

Aa.4 

330.47489 

-59.946312 

(X'12 

(X'lO 

54° 

Aa.5 

330.47529 

-59.946349 

0"18 

0713 

150° 

Aa.6 

330.47546 

-59.946523 

(X'18 

(X'12 

150° 

Ab.l 

330.47571 

-59.943954 

0"11 

0709 

131° 

Ab.2 

330.46741 

-59.947260 

(X'14 

(X'll 

131° 

Ab.3 

330.46852 

-59.944283 

0"11 

0709 

131° 

Ab.4 

330.47515 

-59.946584 

CK'18 

CK'12 

150° 

Ac.l 

330.47489 

-59.943958 

0''25 

0713 

41° 

Ac.2 

330.46669 

-59.947267 

0"25 

0708 

20° 

Ac.3 

330.46912 

-59.943994 

0''30 

0708 

140° 

Ac.4 

330.47441 

-59.946030 

0"12 

0708 

220° 

Ad.l 

330.47537 

-59.943594 

0''28 

0713 

40° 

Ad.2 

330.46685 

-59.946564 

0{'26 

0710 

60° 

Ad.3 

330.46784 

-59.944468 

0''12 

0708 

157° 

Ad.4 

330.47326 

-59.947020 

0{'42 

CK'13 

160° 

Ae.l 

330.47345 

-59.943276 

0''53 

0717 

178° 

Ae.2 

330.46590 

-59.946186 

0"28 

0716 

100° 

Ae.3 

330.46837 

-59.943805 

0''30 

0713 

150° 

Ae.4 

330.47315 

-59.946447 

0"42 

0713 

130° 

Af.l 

330.47417 

-59.943267 

0''52 

0715 

10° 

Af.2 

330.46621 

-59.945961 

0{'39 

0716 

130° 

Af.3 

330.46745 

-59.944289 

0''37 

0713 

123° 

Af.4 

330.47249 

-59.946730 

CK'42 

CK'13 

130° 


line emission can be explained by lensing magnihcation. Ac¬ 
cording to our hducial Lenstool model (see section 4.2.2), 
the magnihcation at the position of the bulge images, and 
variation across the spiral is 1.29lo'i5 for image 1, 0.95lo'o3 
for images 2 and 3, but 1.621 q; 28 for image 4 (with the mag¬ 
nihcation greatest near galaxy N.l). 

The 835.5 nm emission near the galaxies N.4 and N.3 is 
possibly a demagnihed image of bulge Ao and knot Aa. A 
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demagnified image of the bulge is robustly predicted between 
N.2 and N.4, although our fiducial Lenstool model places 
it closer to N.2 at (330.47135, -59.945850). Models allowing a 
demagnified image of Aa stretched towards N.3 are possible, 
but at lower likelihood, and (for reasonable positions of the 
cluster-scale halo) this would be in addition to a demagnified 
image between N.2 and N.4. Since these identifications are 
not robust with current data, and could be merely imperfect 
foreground subtraction near the bright cluster galaxies, we 
exclude them from our strong lensing analysis 

Arc B would be difficult to interpret as a strong lens, as 
previously suggested, because it is at lower redshift than sys¬ 
tem A but greater projected distance from the lens (whereas 
Einstein radius should increase with redshift). Instead, its 
constant velocity gradient suggests that is merely an edge-on 
spiral galaxy, aligned by chance at a tangential angle to the 
cluster, elongated and flexed (our fiducial Lenstool model 
predicts shear 7 = 0.20 ± 0.01) but only singly imaged. The 
two knots of star formation at its north end further enhance 
its visual appearance of curvature. This is consistent with 
the absence of observed counterpart images, and the lack of 
fold structure in its 2D velocity field. We therefore exclude 
arc B from our strong lensing analysis. 


3.3 Modelling the cluster mass distribution 


To model the strong gravitational lens system, we use two 


independent software packages: Grale (Liesenborgs, De Ri- 


jcke & Dejonghe 2006) and Lenstool (Julio et al. 2007). 


The two packages make very different assumptions. Grale 
models the mass distribution using a free-form grid, in 
which the projected density at each pixel is individually con¬ 
strained and individually adjustable; Lenstool interpolates 
a parametric model built from a relatively small number of 
components, each of which has a shape that matches the 
typical shapes of clusters. The two packages also exploit 
slightly different features of the input data. For example, 
both methods match the position of multiply-imaged sys¬ 
tems, but Grale can also match their shape, and use the 
absence of counter images where none are observed; while 
Lenstool can expoit the symmetries of great arcs. 


3.3.1 Grale 


Grale is a free-form, adaptive grid method that uses a ge¬ 
netic algorithm to iteratively refine the mass map solution 


(Liesenborgs, De Rijcke & Dejonghe 2006 Liesenborgs et 


al.|2007 2008a|b 2009 Liesenborgs fc De Rijcke|2012 ). We 
work within a 50^' x 50^' reconstruction region, centered on 

(330.47043, —59.945804). An initial coarse resolution grid is 
populated with a basis set; in this work we use projected 
Plummer density profiles (|h ummer | |1911| |Dejongh^|1987| ). 
A uniform mass sheet covering the whole modelling region 
can also be added to supplement the basis set. As the code 
runs, the more dense regions are resolved with a finer grid, 
with each cell having a Plummer with a proportionate width. 
The code is started with an initial set of trial solutions. 
These solutions, as well as all the later evolved ones are 
evaluated for genetic fitness, and the fit ones are cloned, 
combined and mutated. The resolution is determined by the 
number of Plummers used. The initial coarse resolution grid 


is refined nine times, allowing for more detail in the recon¬ 
struction; the best map is selected based on the fitness mea¬ 
sure. The final map consists of a superposition of a mass 
sheet and many Plummers, typically a few hundred to a 
thousand, each with its own size and weight, determined by 
the genetic algorithm. Note that adopting a specific (Plum¬ 
mer) density profile for our basis set does not at all restrict 
the profile shapes of the mass clumps in the mass maps. 

We use two types of fitness measures in this work, 
(a) Image positions. A successful mass map would lens 
image-plane images of the same source back to the same 
source location and shape. We take into account the posi¬ 
tion, shape, and angular extent of all the images in table 
by representing each image as a collection of points that 
define an area. A mass map has a greater fitness measure 
if the images have a greater fractional area overlap in the 
source plane. This ensures against over-focusing, or over¬ 
magnifying images, which plagued some of the early lens 
reconstruction methods, (b) Null spaee. Regions of the im¬ 
age plane that definitely do not contain any lensed features 
belong to the null space. Each source has its own null space. 
Typically, a null space is all of the image plane, with ‘holes’ 
for the observed images, and suspected counter images, if 
any. The product of these two fitness measures is used to 
select the best map in each reconstruction. 

All lensed images in this cluster arise from extended 
sources (star formation knots within a galaxy). Because of 
that it is hard to identify the center of each image to a pre¬ 
cision comparable to HST resolution. The 073-076 extent 
of Grale points representing some image will be, approxi¬ 
mately, the lower bound on the lens plane rms between the 
observed and predicted images. 

We run twenty mass reconstructions for each image con¬ 
figuration (a limit set by computational time constraints), 
and the maps that are presented here are the averages of 
these. The full range of the recovered maps provides an es¬ 
timate of the statistical error in the mass maps. 


3.3.2 Lenstool 


Lenstool is a parametric method that uses a Markov-Ghain 
Monte Garlo (MGMG) fit with a comparatively smaller num¬ 
ber of mass peaks, but which are free to move and change 
shape. We construct our mass model using one dual Pseudo 
Isothermal Elliptical (dPIE, [Limousin et al. 2005 Elfasdottir 
et al.|200'7 ) halo for the overall cluster, plus a smaller Pseudo 
Isothermal Elliptical halo for each galaxy N.I-4 and N. 6 . 
Each halo is characterised by a position (x, y), velocity dis¬ 
persion (Jv, ellipticity e, and truncation radius Xcut; the clus¬ 
ter halo is also allowed to have a nonzero core radius rcore- 
We set the following priors on the cluster halo: e < 0.75, 
Tcore < 4", and the position has Gaussian probability with 
width g — 2" centered on N.2. For the galaxy halos, we set 
a prior e < 0.45. The position of N.I is of particular scien¬ 
tific interest, and will be well-constrained because it is sur¬ 
rounded by strong lens images, so to avoid any bias, we set 
a prior that is flat within — 5" < x < 3" and — 3" < y < 3" of 
the optical emission. The position of N.2-N.4 will be less well 
constrained, so we set Gaussian prior^with width cr = 075, 


^ We check that the Gaussian priors do not bias the best-fitting 
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centred on their optical emission. The parameters of N.6 
are poorly constrained, because it is faint and far from the 
strong lens systems, so we fix its position to that of its optical 
emission, and fix e = 0 to reduce the search dimensionality. 
The strong lensing data alone provide no constraints on the 
outer regions of the mass distribution, so we manually fix 
^cut = 1000^' for the cluster halo, and rcut = 100^' for the 
galaxy halos, well outside any region of influence. 

As constraints, we use the positions of all multiple im¬ 
ages in table and, following additional symmetries of the 
image, require the za = 1.24 critical curve to pass through 
(330.47113,-59.943529) and (330.46677,-59.945195) per¬ 
pendicular to an angle of 175° and 110° respectively. These 
are indicated in the bottom panel of figure]^ We then opti¬ 
mise the model in the image plane using an MCMC search 
with parameter BayesRate= 0.1 (which allows the poste¬ 
rior to be well explored during burn-in, before converging to 
the best-fit solution; for more details, see |Kneib et al.|1996| 
Smith et al.||200^ [Jauzac et al.|[^14a ) . We assume an er- 


ror of 072 (68% CL) on every position. This choice merely 
rescales the posterior. If we instead assume an error of 0'.'273, 
the model achieves reduced x^/dof=l. 


4 RESULTS 
4.1 Stellar mass 


Spectroscopic redshifts of N.1-N.4 and N.6 confirm that they 
are all 0.099 cluster members. Notably, N.1-N.3 are at 
essentially the same redshift as each other (which is consis¬ 
tent with the mean redshift of all known member galaxies), 
and are projected along a straight line. One the other hand, 
N.4 has a relative line-of-sight velocity ~900km/s. 

In table we interpret the galaxies’ integrated broad¬ 
band fluxes (in individual filters) as stellar masses M* 
via the models of Bruzual & Chariot (2003), assuming a 


Chabrier (2003) initial mass function (IMF), solar metallic- 


ity, and formation redshift Zf = 3. The inferred stellar mass 
of N.l is surprisingly low compared to its stellar velocity 
dispersion, but this could be because it is furthest from the 
cluster core, so its measured flux contains the least contam¬ 
ination from intra-cluster light. 

We also interpret the 3(7 upper limits on the Yla flux as 
star formation rate (SFR), following Kennicutt (1998) rela¬ 
tions converted to a Chabrier (2003) IMF. The red and dead 
ellipticals exhibit effectively zero star formation, consistent 
with their faint UV broad-band fluxes. 


4.2 Total mass distribution 

4-2.1 Grale 

The Grale reconstruction of the projected mass distribu¬ 
tion within the central 30''x 30^' of the 50''x 50^' fitted region 
is shown in the top panel of figure]^ The mean rms offset of 
observed lens images from their predicted positions in the 


position by shifting the priors to the peaks in the posterior and 
rerunning the analysis. Positions all move by less than 0708 in x 
and y, indicating convergence to within statistical errors. 



_J_^^^_ I ^ _ I_L_.^J_I_^ 

0 10 20 
-ARA (arcsec, relative to Nl) 


Figure 4. Top panel: Map of total mass in the cluster core, 
reconstructed using Grale. Green contours show the projected 
mass density, spaced logarithmically by a factor 1.15; the thick 
contour shows convergence n = 1 for Zd = 0.099 and = 1-24 
(Scrit = 1.03 g/cm^). Red dots show local maxima in individual re¬ 
alisations of the mass map. Black dots show cluster ellipticals N.l- 
N.4. Blue circles show the lensed images. Middle panel: Mass after 
subtracting a smooth cluster-scale halo to highlight substructure. 
The thick contour is at = 0. The green (positive) and yel¬ 
low (negative) contours are at = ±0.025, ±0.05, ±0.1, ±0.2.. .. 
Bottom panel: Total mass, as in the top panel but reconstructed 
via Lenstool. The red dashes show the zx = 1.2A critical curve. 
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image plane is (rms^) = 0^/34, but this value is inflated by 
the size of the input images, as explained in Section [5.3.1| 

The mass distribution peaks at (330.47068,-59.945513), 
i.e. (8'.'15, iy73) or r = 8'.'33 from galaxy N.l. The total mass 
of the cluster core inside a cylinder with this radius and cen¬ 
tered on the mass peak is Mc£ = 3.58 x 10^^ M©. The typical 
fractional error in surface mass density in the central region 
is ~ 15%. The total mass projected within l75 = 2.7kpc of 
galaxies N.1-4 is 1.63, 2.02, 1.70 and 1.67 x 10^^ M© respec¬ 
tively. 

To highlight the position of substructures, the red dots 
in the top panel of figure show the local maxima in each of 
the individual, statistically independent realisations of the 
mass map (within the central = 1 contour only), and the 
middle panel of figure shows the mean map after subtract¬ 
ing a smooth, cluster-scale halo. The cluster-scale halo used 
is centered on the peak density, has constant projected den¬ 
sity inside a core of radius 474 = 8.0 kpc, and a projected 
density profile p 2 i:>{r) oc r~^'^ outside this core. 

The most robustly constrained region of interest is near 
galaxy N.l, owing to the proximity of several lensed images. 
Local peaks in individual realisations of the mass map form a 
tight cluster l701±0.39 east-southeast of the peak of optical 
emission, implying that the offset is significant at the 2.6a 
level. The mass within l75 of this location is 1.66 x 10^^ M©. 

Galaxy N.3 is next nearest to gravitationally lensed im¬ 
ages. The closest local mass peak is frequently 1-2" north¬ 
west of the optical emission, and the mass within 1.5" of this 
location is 1.38 x 10^^ M©. Statistically significant structure 
is apparent in the mean mass map, but its presence is less 
robust than for N.l. It may be interesting that galaxies N.l- 
3, the mass peaks closest to N.l and N.3, and the cluster’s 
large-scale diffuse light all lie close to a straight line. Infall 
from preferred directions along filaments is expected {e.g. 


Schaye et al. 2015) — although orbits of galaxies do not 


generally stay radial this far from the virial radius, so it 
may also be coincidence. Galaxy N.4 has no local mass peak 
in many realisations of the reconstruction; the closest set 
of peaks is offset ~ 3" to the southeast. However, the re¬ 
gion around N.4 is not as well constrained by nearby lensed 
images, and appendix [B] shows that the distribution of sur¬ 
rounding mass is sensitive to the assumed identification of 
strong lens images. The position of the mass associated with 
galaxy N.2 is difficult to disentangle from that of the cluster- 
scale halo. 


4-2.2 Lenstool 

Best-fit parameters for the fiducial mass model are listed 
in table and the corresponding mass map is shown in 
the bottom panel of figure (it is uninformative to sub¬ 
tract the cluster scale halo, as its core is flatter than the 
cuspy galaxies, so does not visually affect their position). 
The mean rms offset of observed lens images from their pre¬ 
dicted positions in the image plane is (rms^) = 0726 (41% 
of which is contributed by systems Ae and Af, whose ob¬ 
served position in some multiple images is indeed the least 
certain). With the assumed 072 errors on observed positions, 
the model achieves y^/dof = 49.3/23, and Bayesian evidence 
log^o(^) — —26.4. Some of the Lenstool MGMG samples 
predict fifth and sixth images of knot Ac to the east of Aa.5 


and Aa.6, which are indeed faintly visible as the continua¬ 
tion of image A’s ring. 

The total mass projected within r=8733 of (330.47068, 
—59.945513) is Md = 3.49 x 10^^ M©. The total mass within 
175 of galaxies N.1-4 is 0.65, 1.20, 1.48 and 1.83 x 10“ M© 
respectively. These numbers include a contribution from the 
cluster-scale halo, and are most useful for comparison to 
the Grale results. Integrating the individual components 
of the mass model analytically within a 1.5" radius, exclud¬ 
ing the cluster halo, yields mas ses of 0.72, 0.95, 1.28 and 
1.10 X 10^^ M© (see equation 9 of Limousin et al.|2005 ). The 
mass associated with N.l increases slightly in this calcula¬ 
tion because these measurements are centered on the mass 
peaks. The 68% confidence limits on all of these masses, ob¬ 
tained by propagating the statistical uncertainty (only) in 
the observed image positions through the MGMG sampler, 
are approximately 1%. 

The mass associated with galaxy N.l is offset by 
0776lo‘ i6 fhe peak of its light emissiorj^ The 68% con¬ 
fidence limit, again obtained from the MGMG samples, in¬ 
cludes only statistical uncertainty propagated from that in 
the observed image positions. To estimate the additional 
uncertainty (bias) caused by the relative inflexibility of 
Lenstool’s parametric model to represent a complex mass 
distribution, we retry the optimisation using two alternative 
model configurations. 

Refitting the data using a model with a second cluster- 
scale halo achieves (rms^) = 0725, although with worse 
X^/dof = 48.5/17, and logiQ(L^) = —28.0. In this model, 
the first cluster halo remains between N.2 and N.3, and the 
second cluster halo appears between N.3 and N.4, whose 
masses are reduced by ~ 45%. The mass associated with 
galaxy N.l moves to (-07681°;}^, -0725t°;i^). This 0722 
shift relative to the fiducial reconstruction provides one es¬ 
timate of model-induced bias. 

Second, we can reuse the fiducial model to perform an 
MGMG fit to a mass distribution that has a similar con¬ 
figuration to the real distribution but is not impeded by 
the parametric limitations of Lenstool. We generate four, 
slightly different realisations of mock data by raytracing the 
positions of each lensed image {e.g. A[a-f].l) through the 
fiducial model to the source plane, then back to (multiple 
locations in) the image plane. In this mock data, the true po¬ 
sition of the mass is known, and can be represented perfectly 
by Lenstool. We perform four independent fits, centering 
the priors around the true positions of the mass. The mean 
spurious offset of N.l is 0755 ± 0.11, and the mean spurious 
offsets of N.2-4 are 0755 ±0.09. Gomparing this to the to¬ 
tal ~0776 offsets in the real fit therefore suggests a similar 
~0721 budget for model-induced bias. 

The mean offset of N.2-4 (with the observed con¬ 
straints) is 0762, which could be interpreted to imply a 
characteristic error on N.l of this order. However, these are 
measured with a prior, and their uncertainty is much greater 
than that of N.l because those galaxies happen to be much 
farther from strongly lensed systems — and, in the case of 
N.2, because of degeneracy with the position of the cluster- 


^ This remains as 0781 ^q‘} 2 even if 10% of the mass (accounting 
for the stellar component) is forced to lie on the galaxy; the offset 
is necessary to ensure the observed multiplicity of Aa.5 and Aa.6. 
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Table 3. Parameters of the best-fit, fiducial mass model constructed by Lenstool. Positions are relative to the peak of light emission, 
except for the cluster-scale halo, whose position is relative to the peak of emission from galaxy N.l. Quantities in square brackets are not 
fitted. Errors on other quantities show 68% statistical confidence limits due to uncertainty in the lensed image positions, marginalising 


over uncertainty in all other parameters. 


CTv [km/s] x['] y["] e 0 [°] rcore ["] (rmsQ ["] X^/dof logio(^) 

Fiducial model: 0.26 49.3/23 —26.4 


N.l 

1901®2 

—0 

-0.46l°l° 

0 23+0-15 

101122 

[^0] 

N.2 

2191^1 

— 0 13+0-28 


u.uy_o 09 

174122 

[^0] 

N.3 

254+\l 

u.uy_0 25 

— 0 3C+0-18 

U.C)D_o 29 

0'-^o_o 10 

30l“ 

[^0] 

N.4 

235+34 

—0 99 + 0-30 
u.a»_o 34 

'0-0'J-_o.27 

0 IQ+O-12 
0'-l^-0.09 

121154 

[^0] 

N.6 

islf 

[0] 

[0] 

[0] 

[0] 

[^0] 

Cluster 

620li« 

e-i+fo! 

2.301+®8 

0 70+0-01 

O'- '’0_o.24 

6llt 

30 12+0-23 

o0.iZ_g 43 


scale halo. Contrary to the results from Grale, the mass 
associated with galaxy N.3 is coincident with the position of 
its light emission within measurement error; the mass within 
1/5 of Grale’s offset peak is a lower 1.07 x 10^^ M©. The 
mass associated with N.4 is offset at only marginal statis¬ 
tical significance but, intriguingly, the offset is in the same 
direction as that measured by Grale. 


5 INTERPRETATION 


We have modelled the distribution of mass in the cluster us¬ 
ing two independent approaches: free-form Grale and para¬ 
metric Lenstool. The general agreement between methods 
is remarkable, both in terms of total mass and many details. 

The most striking result is that the mass associated with 
galaxy N.l is offset from its stars, l/01±0.39 east-southeast 
with Grale or O/TGIqJt southeast with Lenstool (lin¬ 
early adding statistical and method-induced errors). That 
the measurements are consistent with each other, and re¬ 
silient to small changes in the strong lens identifications (see 
appendix [b|, supports a robust conclusion that the offset is 
real. To combine the analyses we note that, although they 
start with mostly identical input data, uncertainty on their 
final constraints is dominated by the highly-nonlinear recon¬ 
struction procedures, which are independent. We therefore 
average the best-fit values with equal weight and add their 
errors in quadrature, to infer a combined constraint on the 
offset between mass and light 

5 = (y.'sglo.a? = kpc (68% CL). (1) 


The strong lens conhguration makes galaxy N.l the best 
measured of ah the cluster members, but both Grale and 
Lenstool provide marginal evidence for an offset in galaxy 
N.4, and Grale suggests a similarly unexpected distribution 
of mass near in N.3. 

Interpreting an offset between mass and stars is difficult. 
It could feasibly be caused by different tidal forces or dynam¬ 
ical friction on the different-sized dark matter/stellar haloes; 
partially stripped gas (or unrelated foreground/background 
structures) that contributes to the total mass ( [Eckert et al.| 
|2Q14| [Roediger et al. 2014); or a displacement of the light 
emission due to recent star formation triggered in stripped 
gas. However, to hrst order, tidal forces do not alter the 
peak position. Archival Chandra data also show no sub¬ 
structure near this cluster core. The effectively zero star for¬ 
mation we observe also suggests that broad band emission 


should trace stars that existed before the merger. Instead, 


Williams & Saha (2011) interpreted the offset in terms of 


dark matter’s self-interaction cross-section crDM/'Ri, using a 
toy model of interactions equival ent to an optical depth (see 
their equation (sj^and see also Massey, Kitching & Nagai 


2011 Kahlhoefer et al.|2014 Harvey et al.|2014 ). If stars in 

the infahing galaxy are subject only to gravity, but its dark 
matter also feels an effective drag force, after infall time 
Lnfaii, dark matter lags behind by an offset 


(^(Lnfall) 


GMciMr>M ctdm ,2 

2 2 ^infall? 


( 2 ) 


where Md is the mass of the cluster interior to the in¬ 
fahing galaxy, which has dark matter mass Mdm and cross- 
sectional area at clustercentric radius tdm- Adopt¬ 

ing mean masses from our Grale and Lenstool analyses, 
Mc£ = 3.54 X 10^^, Mdm = 1.19 X 10^^, parameter sdm = 4/1 


following Williams & Saha (2011), and rDM='r = 8. 3, then 


propagating 10% errors on the masses and 0/5 errors on the 
sizes, suggests 


a/m - (1.7 ± 0.7) X 10“ 


Lnfall 

10^ yrs 


cm^/g. 


( 3 ) 


The infah time must be less than 10^° years, the age of 
the Universe at the cluster redshift. Given the lack of ob¬ 
served disruption, cohinearity (and common redshift) of 
N.1-3, they are likely to be infahing on hrst approach from 
a hlament, and moving within the plane of the sky. Thus 
Lnfaii ~ 10^ yrs, the approximate cluster crossing time, and 
assuming this conservative upper bound places a conserva¬ 
tive lower bound on a/m. If any component of the motion 
is along our line of sight, the 3D offset may be larger, so 
our assumption of motion exactly within the plane of the 
sky is also conservative. Using a different set of strong lens 
image assignments (see appendix [b| , we recover the 6 kpc 
offset and correspondingly larger cross-section of jWilliams fcj 
Saha (2011). These image assignments are now ruled out by 


our new IFU spectroscopy, which unambiguously traces the 
morphology of the lens, even through foreground emission 
and point sources in the broad-band imaging. 

We have also measured the mass to light ratios of the 
four central galaxies. Each of them retains an associated 
massive halo. There is no conclusive evidence to suggest 


^ Note that the prefactor in equation (4) of Williams &: Sahli] 
(2011) should be 6.0x10^ rather than 6.0xl0‘x 
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that any of them are more stripped than the others; if any¬ 
thing, the stellar mass of N.l is marginally lower than ex¬ 
pected, compared to both the stellar and dark matter mea¬ 
surements of velocity dispersion. This is the opposite of be¬ 
haviour expected if the dark matter associated with N.l is 
being stripped. 

We have not yet attempted to measure any truncation 
of the galaxy halos, like Natarajan et al. (2009). That mea¬ 
surement will be improved by combining our current mea¬ 
surements of strong lensing with spatially extended measure¬ 
ments of weak lensing and flexion currently in preparation, 
plus multi-object spectroscopic data of member galaxies out¬ 
side the cluster core. 


6 CONCLUSIONS 

We have presented new Hubble Space Telescope imaging 
and Very Large Telescope integrated field spectroscopy of 
galaxy cluster Abell 3827. This is a uniquely interesting sys¬ 
tem for two reasons. First, it contains an unusually corru¬ 
gated light distribution, with four almost equally bright cen¬ 
tral galaxies (not five, as previously thought, because one is 
a foreground star) within the central lOkpc radius. Second, 
a gravitationally lensed image of a complex spiral galaxy 
is fortuitously threaded between the galaxies, allowing the 
distribution of their total mass to be mapped. Because of 
the cluster’s 2 : ~ 0.1 proximity, this can be achieved with 
high spatial precision. We expect this data will be useful be¬ 
yond this first paper, for several investigations of late-time 
dark matter dynamics in the poorly-studied (yet theoreti¬ 
cally contentious) regime of cluster cores. 

We have investigated the possible stripping or deceler¬ 
ation of dark matter associated with the infalling galaxies. 
Most interestingly, combining two independent mass map¬ 
ping algorithms, we find a 1.62^.^]% kpc offset {i.e. 3.3a sig¬ 
nificance) between total mass and luminous mass in the best- 
constrained galaxy, including statistical error and sources of 
systematic error related to the data analysis. Such an offset 
does not exist in isolated field galaxies (or it would have been 
easily detected via strong lensing of quasars). If interpreted 
in terms of a drag force caused by weak self-interactions be¬ 
tween dark matter particles, this suggests a particle cross- 
section a/m ~ (1.7 ±0.7) x 10“^cm^/g. However, the small 
absolute offset < 2 kpc might be caused by astrophysical ef¬ 
fects such as dynamical friction, and it is difficult to conclude 
definitively that real dark matter is behaving differently to 
CDM. Detailed hydrodynamical simulations of galaxy in¬ 
fall, incorporating dark matter physics beyond the standard 
model, are needed to predict its behaviour within a cluster 
environment, and to more accurately interpret high preci¬ 
sion observations. 
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APPENDIX A: COMPARISON OF VIMOS AND 
MUSE SPECTROSCOPY 

Our lEU observations of the same (faint and highly elon¬ 
gated) source with both VLT / VIMOS ( |Le Eevre et al.|2003 ) 
and VLT/MUSE ( [Bacon et al.|[20TQ| provide an early op¬ 
portunity to compare the practical performance of VLT’s 
old and new spectrographs. As is apparent in figure [AT] the 
throughput of MUSE is vastly greater, achieving higher sig¬ 
nal to noise in a much shorter exposure time. Crucially for 
our purposes, the improved flat field stability allows much 
cleaner foreground subtraction, revealing candidates for the 
demagnified fifth image near foreground galaxies N.3 and 
N.4. MUSE^s higher spectral resolution and larger field of 
view are also responsible for two fortuitous discoveries: our 
identification of line emission in lensed galaxy A as the [Oil] 
doublet, which we had previously mistaken for Ha, and the 
identification of galaxy B as an edge-on spiral, which we 
had previously mistaken for a strong gravitational leasing 
arc. Availability of the new MUSE spectrograph has thus 
revolutionised our interpretation of this cluster. 
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Figure Al. False-colour images of the cluster core, manu¬ 
factured from the VLT/VIMOS and VLT/MUSE IFU 3D data 
cubes. The displayed region shows the full field of view of each 
instrument. Red, green and blue channels respectively correspond 
to a narrow-band [Oil] image (created as described in section [3^ , 
Johnson I band and Johnson V band. 


APPENDIX B: ALTERNATIVE LENSED 
IMAGE ASSOCIATIONS 


To draw robust conclusions about any offsets between op¬ 
tical emission and associated dark matter, one needs to ac¬ 
count for all sources of uncertainty in the mass models. 
Statistical uncertainties in the location of the mass peaks 
around galaxy N.l are easily obtained. There are two types 
of additional uncertainties: those arising from the mass re¬ 
construction methods, and those from the image identifica¬ 
tion. We explore the latter in this appendix, by considering 
two alternative image assignments (see figure Bl). 

In the body of this paper, we presented a fiducial model 
using what we believe to be the correct associations - based 
on the surface brightness, colours and sizes of the knots in 
HST imaging, and the extended distribution of [Oil] emis¬ 
sion in IFU spectroscopy. The only questionable region is 
around N.l; the other image assignments are secure (though 
the exact locations of some faint extended images, like Ae 
and Af, are uncertain to within a few tenths of an arcsecond). 
All the image assignment schemes we have tried produce a 
common elongation of the mass distribution along the north¬ 
west to southeast axis, and a statistically significant offset 
between the light and mass of N.l, using both Grale and 
Lenstool. However, the amount and the morphology of the 
offset differs. 

Alternative configuration L represents a minor pertur¬ 
bation of the strong lens image assignments, with Aa.6 
switched to a fifth image of star formation knot Ab.5. In 
the fiducial configuration, Grale reproduces the three im¬ 
ages of knot Aa, but merged into a continuous arc, which 
is not consistent with their appearance on HST images. In 


Table Bl. Changes to the assignments of multiply imaged sys¬ 
tems in table to produce alternative configuration B. Columns 
denote the ID and position of the image, its major and minor 
axis, and the angle of the major axis on the sky, anticlockwise 
from west. 


Name 

RA 

Dec 

Major 

Minor 

Angle 

Aa.4 

330.47447 

-59.946180 

(X'13 

(X'09 

90° 

Aa.5 

— 

— 

— 

— 

— 

Aa.6 

— 

— 

— 

— 

— 

Ab.4 

330.47452 

-59.946347 

0"10 

0"07 

90° 

Ac.4 

330.47440 

-59.946032 

(Y.'or 

(X'06 

70° 



Figure Bl. Two alternative configurations of lensed image asso¬ 
ciations near cluster galaxy N.l. There is a surprising high density 
of point sources, so it is not obvious which sources correspond to 
which lensed images. Configuration L is a small perturbation of 
the fiducial model in the main paper. Configuration B is based 
upon Lentool’s best-fit mass model if no identifications are ini¬ 
tially made in this region. However, it is at odds with narrow-band 
imaging of the source, created from line emission in our VLT IFU 
data (see figure]^. 


configuration L, Grale fin ds a mass 4.44 x 10^^ M© within 
10'' of the peak (see figure B2), which is better able to re¬ 
produce the observed distribution of luminous sources near 
N.l. However, it introduces an extra mass clump ~ 6" south 
of N.l that has no identifiable optical counterpart (see fig¬ 
ure B3). With this configuration, Lenstool produces in¬ 


consistent predictions of the multiplicity of knot a, and an 
overall fit marginally worse than the fiducial model. In the 
best-fit model, the total mass is 4.48 x 10^^ M© within the 
same aperture as above. The positions of the mass peaks 
move only within statistical errors compared to the fiducial 
model (see table B2), for this and similar minor perturba¬ 
tions to the image assignments. 

Alternative configuration B had been our preferred con¬ 
figuration before we obtained any VLT integrated field spec¬ 
troscopy. If star formation knots Aa.4-6, Ab.4 and Ac.4 are 
left unassigned, and the fit is optimised without any con¬ 
straints in that region, the best-fit model converges to an 
assignment like that in table |B1[ reflecting a morphology 
illustrated in figure |B1| The images immediately south of 
N.l in the broad-band images are unaccounted for as strong 
lens sources, but point sources are frequent in this region, 
so they could be a chance superposition. Using these assign¬ 
ments, Grale finds a total mass 6.06 x 10^^ M© within 10" 
of the peak. The mass near N.l becomes a Tail’ extending 
to the southeast (see figure B4). This is particularly inter- 
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Table B2. Parameters of the best-fit mass model constructed by Lenstool, using alternative associations between lensed images. 
Positions are relative to the peak of light emission, except for the cluster-scale halo, whose position is relative to the peak of emission 
from galaxy N.l. Errors show 68% statistical confidence limits due to uncertainty in the lensed image positions. 
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Figure B2. Total mass profile of the cluster core as a function of 
projected distance from the peak of the Grale mass map, for the 
fiducial model versus the alternate configurations of lensed image 
associations. Circles along the bottom show the radial positions 
of cluster member galaxies N.1-4, and vertical lines mark the 
positions of lensed images. 


results is reassuring: the star formation knots near N.l were 
not resolved in the ground-based imaging. However, our IFU 
data shows this image configuration to be incorrect. [Oil] 
line emission is observed immediately south of N.l in fig¬ 
ure inconsistent with the configuration B morphology il¬ 
lustrated in figure |B1| Note that the observed line emis¬ 
sion flux has a large spatial gradient, and is detectable even 
through the obscuration of N.l because (according to our 
fiducial model) the gravitational leasing magnification is 
53% greater at the position of knot Aa.4 than at the posi¬ 
tion of knot Ad.4. The change in conclusion from |Williams| 


& Saha (2011) to this work thus emerges equally from both 


our new high resolution HST imaging plus our new VLT in¬ 
tegrated field spectroscopy, and demonstrates the discrimi¬ 
natory power of the combined observations. 


This paper has been typeset from a TgX/ HTg^ file prepared 
by the author. 


esting because in some models of particle physics (Kahlhoe- 
fer et al.|2014 ), non-zero self-interaction cross-section would 


indeed disperse the dark matter, rather than simply offset¬ 
ting it from the light. A similar extension would also prob¬ 
ably be expected in the case of dynamical fraction. Using 
configuration B, the best-fit Lenstool model with mass 
4.48 X 10^^ Mq is apparently better than the fiducial model, 
with (rmsi) = 0.23'h 

Most notably, configuration B’s 3'.'351 q; 82 offset be¬ 
tween the mass and light of N.l reproduces the result of 


Williams & Saha (2011). This recovery of the ground-based 
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Figure B3. Map of total mass in the cluster core, as in figure^ Figure B4. Map of total mass in the cluster core, as in figure^ 

but using Alternate Configuration L to associate lensed images. but using Alternate Configuration B to associate lensed images. 


© 2014 RAS, MNRAS 000, [iHT^ 



















































